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AN INCREMENTAL PROCEDURE FOR SOLUTION OF NONLINEAR 
PROBLEMS WITH APPLICATIONS TO PLATES AND SHELLS 


By Manuel Stein 
Langley Research Center 

SUMMARY 

An incremental step-by-step procedure is described for obtaining solutions to 
nonlinear boundary value problems. This procedure is especially suited for computer 
solution of problems cast in numerical form; it is used in this paper to obtain numerical 
solutions for the clamped plate with uniform lateral load, for the stretched simply sup- 
ported plate with a central concentrated load, and for a spherical cap with internal pres- 
sure and an outward acting central concentrated load. There is good agreement between 
the present results and exact results which have been obtained previously. 

INTRODUCTION 

Linear equations are often inadequate for use in the analysis of space vehicle and 
aircraft plate and shell structures. Some design problems involve lateral deflections 
many times the plate or shell thickness such that the neutral surface strains depend non- 
linearly on lateral deflection. Other problems may involve nonlinear material properties 
such as those due to temperature rise or to plasticity. The design of such structures may 
be concerned primarily with strength or with deformation. Useful estimates of the 
deformations can usually be found from rough approximate solutions of the nonlinear 
equations. To obtain meaningful stress distributions, however, much more accurate 
solutions of the nonlinear equations are required. 

Many methods are available for solving nonlinear problems. (See ref. 1.) Analyti- 
cal solutions of the differential equations of the problem are sometimes possible through 
special devices or inverse procedures. Direct methods include iteration, Newton T s 
method, and methods which make use of trial and error. For problems involving a small 
parameter, a regular perturbation method may be employed; for problems involving a 
boundary layer, methods of asymptotic integration may be employed. Weighted residual 
techniques are also available which include the energy method, the Galerkin method, col- 
location, and the method of least squares. These available methods (for a fairly complete 
survey with examples, see ref. 1) have certain limitations and often are applicable to 
rather restricted types of problems. 



The purpose of this paper is to present a procedure which would have relatively 
general applicability to boundary value problems. In particular, the idea of incrementing 
a control parameter is used in such a way that a solution may be obtained for a broad 
spectrum of practical values for this parameter. 

By using this procedure, solutions are obtained for the deformations and stresses 
in the following axisymmetric problems: a clamped circular plate under uniform load 
(sketch 1), a simply supported circular plate with a central concentrated load and uniform 
edge forces (sketch 2), and a simply supported spherical cap with central concentrated 
load and with internal pressure and corresponding edge forces (sketch 3). The results 
obtained for these problems are in good agreement with "exact" theoretical results and 
with the end results given by membrane theory. Although all the examples in the present 
paper are for axisymmetric problems, the method basically is not limited to such prob- 
lems and can also be applied to problems of asymmetric loading and deformation. 



Sketch 1 



SYMBOLS 


a radius of plate or chord radius of shell 

D plate stiffness, 

D nondimensional 

E Young's modulus of material 

a 2 

H shell rise, — 

’ 2R 


Et J 




12 


tiffness parameter, 


d(l 


M 2 ) 

Eta 2 Q 2 / 3 


2 


Mi 


radial plate bending moment 


q 

Q 


number of stations considered 
concentrated lateral load 
uniform lateral load 

P(l - M 2 ) 


load parameter, 


7rEta 


or 


qa 


a - « z ) 


Et 


radial distance 


R 

t 

T 

u 


radius of spherical shell 


thickness of plate or shell 


^ral 1 " M 2 ) qR(i - u 2 ) 

nondimensional edge radial stress, -j - — or — * — y 

EQ 2 / 3 2EtQ 2 / 3 

radial displacement 


w 


deflection 


w f 


central deflection 


function defining the spherical shell considered 


0 

X 


slope, w' 

curvature parameter, X^ = ^ 2a 

rV 

Poisson's ratio of material 


radial stress 


'ro 


-’ra 


central radial stress 


qD' 

edge radial stress [for spherical shell, a ra = 


Primes indicate differentiation with respect to independent variable. 
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DESCRIPTION OF INCREMENTAL PROCEDURE 


The procedure makes use of a one-term Galerkin solution to obtain approximate 
results and a subsequent linearization to improve these approximate results. 

The problems examined herein can be written in the general form 

AL(w) + N(w) = l (1) 

where L is a linear differential operator, N is a nonlinear operator, A is the con- 
trol parameter (or a measure of several parameters but not necessarily large or small), 
and l is a loading term. The choice of the coordinate system and the details of the 
method of solution depend on the particular problem solved; here, let w = w(x). It is 
presumed that a solution is available at each step which satisfies this equation and the 
boundary conditions for a given value of A, for example, A^ n "^. Call this solution 
f(n-l). For a neighboring value of A = A( n ), let w be approximated by 


w 


(n) 

approx 


c (n) f (n-l) 


( 2 ) 


Using equation (2) as a one-term expansion function in the approximate Galerkin solution 
of equation ( 1 ) leads to the following equation which determines the constant c^ when 
f(n-l) 

is known: 

J f( n_ 1 )[c^A^L(f^ n “ 1 ^) + N(c^f^ n_1 ^)] dx = J f( n_1 ^ dx (3) 


(See refs. 2 and 3 for a discussion of the Galerkin method.) Note that equation (3), with 
limits on the integrals specifying the boundaries of interest, is a nonlinear algebraic 
equation in c^ (of same degree as the differential eq. (1)). It is important that the 
proper root of the nonlinear algebraic equation for c^ is chosen as dictated by the 
physical problem. 


Now set 


With F^ assumed to be sufficiently small compared with i vl ~ x> , a linearized form of 
equation ( 1 ) may be written as 

c (n) A (n) L ^ F (n)^ + ^ F (n), f (n-l)^ c (n)^ = l _ c (n) A ( n ) L ( f (n-l)j _ N ( c ( n ) f (n-l)j ( 5 ) 

where N is the part of Njc^(f^ n- '*'^ + F^) which contains F^ linearly. Equa- 
tion ( 5 ) and the boundary conditions determine F^ when c^ and f( n- l) are known. 
With equation (4) defining f^, the procedure can be repeated for the next value of 
A = A^ n+1 \ The present method uses equations (2) to (5) to advance from a solution 


w (n) = c (n)^ f (n-l) + F (n)j = c (n) f (n) 


( 4 ) 


(n- 1 ) 


available at one value of A to a solution at any desired value of A. For each value of 
A - that is, for each step - a check must be made to insure that the linearization assump 
tion is fulfilled (this validates the solution) and, if not, the solution for this step must be 
discarded and the small increments in A must be reduced in magnitude. Thus, the pres 
ent procedure combines the incrementing of the control parameter A with the one-term 
Galerkin method to effect accurate solution over the complete range of A. At each step 
the solution is validated a posteriori. 

For problems of the form of equation (1) the starting point might be the solution of 
the linear equation obtained from equation (1) by neglecting the nonlinear term. Call this 
solution f^. Then, for a large value of A = A^) (an appropriate a(D may be found 
by trial and error) the procedure begins with equations (2), (3), (4), and (5) in which n = 1 
for this first step. In the next and all ensuing steps, A is decreased in small incre- 
ments to A( 2 ), to a< 3 >, and so on to A©, where A© is the A of interest, or until 
A approaches zero, if it is desired to cover the complete range. The procedure just 
described is especially suited for high-speed digital computer solution of problems cast 
in numerical form. 


APPLICATIONS 

In order to illustrate the use of the incremental procedure, solutions are presented 
for two problems. The first problem is the determination of the stresses and deforma- 
tions of a clamped circular plate under uniform lateral loading. This problem was 
chosen because "exact" solutions existed for part of the nonlinear range and for the mem- 
brane case. The second problem considers a simply supported spherical cap subjected 
to internal pressure with corresponding edge forces and a central concentrated outward 
acting lateral load. For zero curvature this problem reduces to that of the flat plate with 
in-plane tension and a concentrated central load. Existing solutions are available for this 
problem only for the end point of the flat membrane. 

In the first problem the critical behavior is near the edge, whereas in the second 
problem the critical behavior is near the center. Both problems were solved numeri- 
cally on the digital computer by replacing derivatives by differences and integrals by 
sums. The conversion of the equations into numerical form is described in the appendix. 


Clamped Circular Plate With Uniform Lateral Load 

Sketch 1 illustrates the configuration and the loading considered. The large deflec- 
tion (nonlinear) differential equations to be solved (see, for example, ref. 4 or 5) are 


D (s" + ? ? ' - 7 ») - 7^2 4' + * “ + i = ' f 

u " + F u '-^" + g ( S ' +1 ^ £ ^)' 0 


( 6 ) 


I 





iniami 


where the dependent variables are the slope 0 = w' and the radial displacement u and 
the primes indicate differentiation with respect to the independent variable r. 

The boundary and continuity conditions considered are 

0(0) = 0(a) = u(0) = u(a) = 0 (7) 

and, if it is required to find the deflection w, there is the additional condition w(a) = 0. 

The numerical equivalent of the method of solution described in the previous sec- 
tion is to be used. Since the region of an expected rapid change of slope is in the neigh- 

r 2 

borhood of the edge, the transformation x = is made to permit equally spaced stations 

a^ 

in the numerical procedure but more closely spaced stations near the edge in the physical 
problem. It is also convenient (but not necessary) to introduce new dependent variables 
B and U according to 





With these changes equations (6) become 

4DxB" - Eta ~ B(2U' - 1 ^ U + i-B 2 ) 

1 - fx2 \ x 2x j 

4xU" + b(2B' - B) = 0 



J 


( 8 ) 


(9) 


f 

\ 


where the primes now indicate differentiation with respect to x, and the boundary condi- 
tions in terms of the new variables are 

B(0) = B(l) = U(0) = U(l) = 0 

From the second of equations (9) and the boundary conditions of U, 

U =±—JL x f 1 ^B 2 dx +i-±Ji(x f 1 i B 2 dx - f * i B 2 dx) (10) 

8 J x x 2 8 ^ J x x J Q x y 

so that the first of equations (9) becomes 



(ID 


P 
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I 


1 


Let an approximate solution be given by 

B = cf (12) 

where f is a known function of x which satisfies the boundary conditions and c is a 
constant that is to be determined by the Galerkin equation. The differential equation (11), 
after substitution for B from equation (12), becomes 


1 ~ U?. ( f 1 ~ f 2 dx + A f * A f 2 dx) + t 1 - f 1 A f 2 dx 
8 Wx x Jq x J 8 J 0 x 


4Dcxf" - c f 
1 - M 2 


To get equation (11) in dimensionless form, divide by qa ; thus, 


qa°x 


(13) 


4Dcxf" - c 3 f 


1 - u 2 




X 

2 


(14) 


where 


D- D(l- M 2 ) t 2 

Eta 2 Q 2 / 3 12a 2 Q 2 / 3 


5 = Ql/3 

_ qa(l - M 2 ) 

Et 


The following Galerkin equation is used to determine c: 



(15) 

This equation corresponds to equation (3) if A = D. Now with c known let 

B = c(f + F) (16) 


with F a function assumed to be sufficiently small, so that squares and higher powers 
of F can be neglected, compared with a known function f. Consequently, a linear dif- 
ferential equation may be written for F: 
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Equation (17) corresponds to equation (5). The numerical procedure used on the digital 
computer to solve these problems is described in the appendix. 

Results for deflections and midplane stresses at the center of the plate are plotted 

as a function of D = ~ in figures 1 and 2. There is very good agreement 

Eta^Q 2 /^ 

between the present numerical results and the exact solution by power series of Way 
(ref. 6) using the same equations as used in the present paper. The present results 
extend the solution of Way to thinner plates or higher loads (smaller D) and, in the limit 
as D — 0, approach the exact membrane solution of Hencky (ref. 7). 

For comparison, the one-term Galerkin solution of the nonlinear equilibrium equa- 
tions using the function obtained from linear theory was determined; it gives only a fair 
approximation for the complete range of D as shown in figures 1 and 2. 




EtaV 3 


Figure 1.- Deflection at the center Figure 2.- Radial midplane stress at the 

of a uniformly loaded clamped center of a uniformly loaded clamped 

circular plate. Q = & circular plate. Q = 

In figures 3 and 4 radial distributions of deflection and in-plane stresses are plotted 
for three different values of D. The deflection distribution for the highest value of D 
as obtained from linear plate theory is also plotted. These results indicate that the 
deflections according to linear theory are everywhere larger than those given by nonlinear 
theory. No unusual changes in the shapes of the curves occur as D changes. Of course, 
in the limit as D approaches zero — that is, for the pure membrane - the boundary 
conditions change from clamped to supported. 
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a 


Figure 3.- Deflection distribution of a uniformly loaded 

03 (l ~ m2) 

clamped circular plate. Q = M N ^ p \ 


Figure 4.- Radial midplane stress distribution of a uniformly 

_ qa(l - p 2 ) 

Et 


loaded clamped circular plate. Q 


Shallow Spherical Shell 


The problem of a simply supported spherical cap with central concentrated load and 
with internal pressure and corresponding edge forces is illustrated in sketch 3. For zero 
curvature the problem reduces to that of a simply supported circular plate with a central 
concentrated load and uniform edge forces (sketch 2). The large deflection (nonlinear) 
differential equations to be solved (see ref. 5, for example) are 




D (e" + ? ? - ji *) ' + Z °)(“' * r u + & + 5 <> 2 ) = f + 5F 


u " + 7 u ' " 72 u + z ° /3 ’ + 0 Z 3 ' + z o + + 2z o) = 0 


j 


(18) 


where 


H 


being the shell rise 




R the radius of curvature, and again the dependent 


variables are the slope /3 = w' and the radial displacement u. Primes indicate differ- 
entiation with respect to the independent variable r. The boundary and continuity con- 
ditions are 
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m = u(o) = o 
4' + r f),« - 0 


E (u'+fu+z^ + V) -Oj 


1 - ll“ 


r=a 


-ra 


( 19 ) 


with the additional condition on deflection at the edge 

w(a) = 0 

For convenience a new dependent variable U is introduced to get homogeneous 
boundary conditions: 


u _ a ra^ ~ ^ £ , tj 
a E a 


(20) 


The region of rapid change of slope for this problem is expected to be near the center; 
therefore, it is expedient to change the independent variable to y where 


y 



With these changes equations (18) become 




u ra (l - n 2 ) ! 


E 


y + f u 




= ^l y 4 + Pa 

2 y + 2 tt 


( 21 ) 


4y 


tU" +4u’ - - § yV + y/3 




= o 

J 


where the primes now indicate differentiation with respect to y. Poisson’s ratio ju is 
set equal to one-third and both sides of equations (21) are multiplied through by y to 


give 



-OrataV^-^y 2 ) 

^ y2 )fcW y2/3D) '-^^ + ^ 2 ' 


qa 3 c Pa 

— y * w y 


( 22 ) 


i(yu’)' - i U - S y VS(,10/3pJ + 1 y 2/S( y «/Sp2)’ . „ 
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qD 

If the edge stress is identified as the stress due to the pressure forces, then a ra = 

and since 3 = JL there is a cancellation with the q term on the right-hand side of the 
2R 

first of equations (22). The boundary conditions in terms of the new variable are 


0(0) = U(0) = 0 

(y 2 / 3 f)y=l - 0 




(23) 


= 0 


y=i 


j 


Because of the different transformation used in this application compared with that used 
in the clamped plate application, it is not convenient to integrate the second of equa- 
tions (22) directly. With f and g known functions of y and c^ and C2 constants 
that are to be determined later by the Galerkin equations, let 


0 = c i f 
U = c 2 g] 

Upon substitution of equations (24) into equations (22), the following equations result: 


(24) 




f - Cl y 3 l 




Eta 

r 

1 - ^ 


L f _ a 2\ 

c 2 f 2/3 

V a 3. y 2,2 

l 0 ! 1 R 7 / 

2y^73\ 

g; r c l y 1 + 2 1 


= 

2tt < 


(25) 


c ,g(W)' - ± g ] - i c iy V3( y 10/3 f )’ . 2|! y V3( y 4/ 3f 2j . 0 
To get in dimensionless form, divide by Pa/V; thus, 

D5ij|(y f ')’ -if]- Tc^f - (c x f - Vy2)j^ y 4 / 3 (y 2 / 3 g) 

-D V VVf + ^yV 




y 

2 


.gW - i g ] - - !|! y V3( y 4/ ¥ )' . 


< 26) 
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where 


D = ~ ^ 2 ) = t 2 

Eta 2 Q 2 / 3 12a 2 Q 2 / 3 

T _ qR(l - /I 2 ) _ CT ra(i ~ M 2 ) 
2EtQ 2 / 3 EQ 2 / 3 


4 QJ 

X 4 = 12 -4- 


Q - P(l - M 2 ) 

Q 71 Eta 


Cl = Ql/3 

- _ c 2 

C2 ‘ q2/3 


The following Galerkin equations may be used to determine Cj and c 2 : 


£ f fogw - if] - Tc^y^f . (c,, - D 1/2 x 2 y2)^ y V3( y 2/3 g )' . 5 1/2 *V £ 


- T- y3f2 / dy " I -C £y dy 

£ g L|(y gr -»,]-> D^V'V 073 *)' - y^/^jL . „ 


Equations (27) are treated in a manner similar to that for equation (3) with A = D. 

Let 

(2( 

U = c 2 (g + G)J 

with functions F and G assumed to be small on the average compared with functions 
f and g, respectively, and with Cj and c 2 known constants. With the aforemen- 
tioned assumptions, the following linear differential equations can be written for F 
and G: 


12 



DcjgfrF’)’ - | F] + |D 1/2 X 2 5 2 y 10 / 3 (y 2/3 G)' - Tc/p - Dc^Vf + SD^VcjVfF 

*¥ y 4/3 (y 2/3 g)’ . | 


= \ y - DCjliyf')' - i f] + Tc lY 3 f + (cjf - D 1 / 2 X 2 y 2 ) 

|y 4 /3( y 2/3 e )’. B l/ V - ilr5f .!i!y3^ 

c 2 [i(yG')' - i G] - J D 1/2 X 2 c iy 2 / 3 (y 10 /3F)' - y 2 / 3 (y 4 /3 fF )’ 


= -c 2 |(yg*)’-ig] + |i 

b 1/ \h l yVyVh) . ,V3( y 4/V)' 


Equations (29) correspond to equation (5). 

The numerical procedure used on the digital computer to solve this problem is 
described in the appendix. Results for deflections and the midplane stresses at the cen- 

— D (i _ u 2) cr ra (l — 

ter are plotted as a function of D = — *-5 — £7^ for values of T = s-70 — between 

Eta 2 Q 2 / EQ 2/d 

0 and 1 for a flat plate in figures 5 and 6. The center deflection for T = 0 and T = 1 
for a flat plate is compared with corresponding results for spherical caps in figure 7. 
Similar comparisons for center midplane stress are shown in figure 8. Similar com- 
parisons of deflection distribution and radial stress distribution are shown in figures 9, 
10, 11, and 12. 



Pfl - M 2 ) 

r,.2«2/3 


Figure 5.- Deflection under a concentrated lateral 
load of a simply supported circular plate with 

PC 1 _ m2) 

edge stress o ra . Q = — -rf-, - 



0(1 ■ u 2 ) 

r,2 n 2« 


Figure 6.- Radial midplane stress under a concentrated 
lateral load of a simply supported circular plate with 

edge stress o ra . Q = • 
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The results show the trend 
toward membrane behavior as the 
thickness decreases or the load 
increases (as D — 0). In figure 8 
the flat-plate deflections for 
T = 0.6, 0.8, and 1.0 approach the 
nonlinear membrane values 
obtained from Schwerin's results 
as D — 0 (ref. 8). For values of 
T lower than 0.5, membrane 
theory gives infinite deflections; 
the stress at the edge is not enough 
to prevent this unlimited deflection 
(the membrane is not in static equi- 
librium). Also, in figure 10 the 
deflection distribution curves for 
D = 0.0013 (and T = 1) illustrate 
the present results when mem- 
brane behavior predominates for 
both the flat plate and the spherical 
cap. 

Figures 11 and 12 show that, 
although the flat plate and spherical 
cap have roughly the same mid- 
plane radial stresses under the con- 
centrated load, away from the con- 
centrated load the radial stress in 
the flat plate decreases and the 
radial stress in the spherical cap 
first increases appreciably and then 
decreases. 



Figure 11.- Radial midplane stress distribution for a circular plate and a 

°ra(l - |j 2 ) 

spherical cap with zero edge stress. 2/3 


0(l - M 2 ) 



Figure 12.- Radial midplane stress distribution for a circular plate and a 
spherical cap with prescribed edge stress. ° ra ^ 2 / 3 ^ ^ = t 


CONCLUDING REMARKS 

A procedure for the solution of nonlinear problems has been presented which 
involves step-by-step incrementing of a control parameter in such away that solutions 
may be obtained for a wide spectrum of practical values of the parameter. The proce- 
dure has been applied to the solution of several axisymmetric problems in plates and 
shells having moderately large deformations. In these problems the control parameter 


15 



is the flexural stiffness. Results are obtained for the complete spectrum of parameter 
values, from the linear plate or shell solution to the nonlinear membrane solution. These 
solutions provide an evaluation of the suggested procedure in three ways. First, a 
posteriori examination of the nonlinear terms at each step of the procedure shows that 
they are sufficiently small to justify linearization and indicates the appropriate size of 
increments in the control parameter. Second, close agreement is obtained with previ- 
ously known linear and nonlinear plate and shell solutions. Third, the plate and shell 
solutions approach closely known nonlinear membrane results for the limiting case of 
the flexural stiffness approaching zero. 

It is believed that this procedure is fast, accurate, more flexible, and more widely 
applicable (with appropriate modification) than many of the methods of solution of non- 
linear boundary value problems now in existence. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 11, 1966, 

124-08-06-01-23. 
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APPENDIX 


GOVERNING EQUATIONS IN NUMERICAL FORM 

The equations of the section entitled "Applications” are herein cast in numerical 
form in order to proceed with their solution on the digital computer. Stations are meas- 
ured out radially from the center of the plate or shell in terms of x or y, and these 
stations are taken to be equally spaced and are numbered from s = 0 at the center to 
s = S at the edge. The following (central) difference approximation for the first deriva- 
tive is used: 



where s indicates the station and e the distance between stations. With the differen- 
tial relations replaced by such differences, the differential equations are then defined for 
each station. Integrals are approximated by finite sums according to the trapezoidal 
rule. Of course, the equations written for the stations near the center and near the edge 
must be modified to account for continuity and boundary conditions. Thus the linear dif- 
ferential equations are replaced by a set of simultaneous linear algebraic equations in 
terms of the unknown displacements at the various stations. 

Clamped Circular Plate With Uniform Lateral Load 

Equations (15) and (17) are now expressed in numerical form. The algebraic rela- 
tion determining c corresponding to equation (15) can be written in matrix form as 

ffj(4DS[x][d] - S S [l][M])|f| = i[fj|x| (Al) 

where 
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Mj 



M 2 



and vertical lines 
matrix — that is, 


f 


1 




indicate column vectors and the half brackets indicate a row 


*1 

*2 


f S-l 


Also, 


and 


where 


1 

2 


x 


1 

S 


S-l 


\| - |£l f 2 • * ' f S-l) 


M 




1/2 1 



1 

1/2 1 
1/2 
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ffl 



. . 1 

. . 1! 


i 


i 

X 


= S 


1/2 


1 

S-l 


1 



1/4 


1 

(S-l) 2 


and the superscript T indicates the transpose. 

The set of simultaneous algebraic equations replacing equation (17) can be written 
in the matrix form as 

- c 3 ([M] ♦ [l][o]jJ |F| - - i|x| - (45c[x][d] - c 3 0[m]) |«| (A2) 

To obtain a solution for a particular case, it is necessary to specify Poisson's 
ratio /j., number of stations S, and initial values of |f|, c, and D. Results are pre- 
sented for / u = 0.3 and S = 12. (It was found that accurate results were obtained for 
S = 12.) The initial values of | f | were obtained from the linear (continuous) solution 
so that 

f s = § - p (s = 0, 1, 2, . . ., S) 

Trial and error suggested a starting value of D of 0.056; then from a Galerkin solution 
using this linear expression, c = 1.032. Insert these initial values in equation (A2) and 
solve for |f|. If the elements of | f| are small compared with | f|, this may be taken 
to be a valid solution so that all the quantities of interest are determined at this value of 
D. For this problem the criterion used to satisfy the linearization assumption was to 
require that the sum of the magnitudes of the components of | F | was less than 4 per- 
cent of the corresponding sum of those of | £ j. From a study of the Galerkin solution a 
change in D of 0.005 was chosen so that, with the new value of D (D = 0.051) and |f| 
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(i.e., the sum of the old |f | and the |f| as just found), a new value of c is deter- 
mined from equation (Al). With these new values of c, D, and |f|, anew |f| is 
determined from equation (A2). This procedure is repeated until results for c, jf|, 
and |F| are obtained for the D value of interest. 

When each new |F) is determined, it should be checked to see whether it is small 
compared with the |f| used to obtain it. If it is not small, the change in D should be 
decreased from the previous step and this step should be recalculated. 

The quantities of interest which may be desired at each value of D may be cal- 
culated from the numerical equivalents of equations (8), (10), and (12) and from /3 = w’. 
In particular, the deflection and the midplane stresses at the center of the clamped cir- 
cular plate were calculated by using Simpson’s rule for integration. The center deflec- 
tion is 


w 0 


aQV3 2S 


1 

x 


f 


(A3) 


The radial midplane stress is given by 


ffr= r^( u ’ + £ u+ H 


so that the stress at the center is 


°ro( ] - c 2 il - u 2 IViri 
EQ 2 /3 " S \ 8 




2 r . i Pi 


(A4) 


where 


T S J = f[ 6 ’ 1 4 > 2> 4 > 2 > • • - 4 > 2 > 4 ] 

[s j = f[ 4 > 2, 4, 2, . . ., 4, 2, 4 ] 

Center deflection and midplane stress as given by equations (A3) and (A4) are plotted as a 
function of D in figures 1 and 2. Figures 3 and 4 present deflection and stress distribu- 
tions obtained for particular values of D. 
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Shallow Spherical Shell 


Equations (27) and (29) are now expressed in numerical form. The algebraic rela- 
tions determining c^ and Cg corresponding to equations (27) can be written in matrix 
form as 
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1 

2 



L S J 

and the M’s required for the diagonal matrix appearing in fA^jl are given by 



The set of simultaneous linear algebraic equations replacing equations (29) is 



where 



and 
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and 



L F sJ 

^ 1/2 

G 3/2 



To obtain a solution for a particular case, it is necessary to specify the number of 
stations S, the tension at the edge T corresponding to the pressure, the curvature 
parameter X, and an initial value of D. The initial value of D chosen should be in the 
range where linear theory is expected to be valid. The initial values of f and g are 
found from the following relations: 





(A7) 


0 


J 


less the last three terms (the terms nonlinear in c^ and Cgj. 

first term in [a^^], and |a. 2 jJ is only the first term in jA^jj . Thus, 

the initial values of | f | and |g| are the numerical solutions to basic equations (26) 
with nonlinear terms omitted. Use these initial values with the initial value of D in 
equations (A5) and solve for Cj and Cg. 

Equations (A6) may be solved for | F | and | G | for these initial values of c«, Cg, 
D, |f|, and j g|. If elements in the |F | and |G | are small compared with the largest 
elements in |f| and jgj, then the solution for the initial value of D is complete. If 


where {Aji] is 
^ 12 ] is only the 
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not, start again in equations (A5) with a larger value of D until the solution for the ini- 
tial D is complete. Then smaller values of D can be considered. 

A new (smaller) D for the next step is obtained by subtracting a prescribed 
increment from the previous D. In equations (A5) use the new value of D and the new 
values for |f| and |g| equal to the previous value of |f| plus |F| and the previous 
value of J g | plus |G|, respectively, and solve for new Cj and Cg. To obtain c^ 
and Cg, it will be necessary to solve a cubic equation and to determine the proper root. 
Use the new values of D, Cj, and Cg and the new jf | and | g| in equations (A6) and 
solve for new |F| and |G|. If the elements of | F | and | G j are of permissible size, 
the solution is now complete for the new D. With each incremental decrease in D, this 
cycle is repeated until a solution is obtained for the D desired. 

The quantities of interest which may be desired at each value of D may be cal- 
culated as in the clamped plate problem. In particular, the deflection and the midplane 
stresses at the center of the shallow spherical shell were calculated by using Simpson f s 
rule for integration. The center deflection is, therefore, 


w 0 ^ c ln -i 

^73 = ii^’ 2 ’ 4 ’ 2 ’" ” 4 ’ 2 ’ 4 ’ 4 ] 


aQ 


The equation determining the midplane radial stress is 


*1 

2f n 


Si c 


(A8) 


+ *o'e 4 e 2 ) 


so that at the center a r is given by 

CT r o(l - M 2 ) - s 2yf3\ 2 / 3 _ 1 

e q2/3 2 jU) 3/2 \2J 


g l/2 2(2) 5 / 3 


( 5\ 2 / 3 / 

3\ 2 /3 1) 

(2) g 5/2-(' 

2) g 3/2 L 


(A9) 


In figures 5 and 6, these quantities are plotted as a function of D for values of T 
(0 s T 5 1) for a flat plate. The center deflections for T = 0 and T = 1 for a flat plate 
are compared with corresponding results for a spherical cap in figure 7. Similar com- 
parisons for center midplane stress are shown in figure 8. Comparisons of the deflection 
distribution for T = 0 and T = 1 for a flat plate and the corresponding spherical shell 
are shown in figures 9 and 10, respectively. Similar comparisons of the radial stress 
distribution are given in figures 11 and 12. 
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